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Introduction 

The phenomenon of vortex breakdown is of considerable importance for flows under 
critical conditions about delta wings. The leading edge vortices on the leeward side 
of the wing can be termed slender until breakdown is about to occur. The 
slenderness of the vortices is characterized by the ratio of the core radius R and 
the length L of the undisturbed vortex, i.e. R/L«l. As long as the slenderness 
condition is satisfied, the equations of motion can be simplified considerably for 
the description of the flow in the vortex. The integration of the simplified 
equations is then possible for large Reynolds numbers. The breakdown itself must 
be analysed by solving the full Navier-Stokes equation. With the slender vortex 
approximation the overall computation time can be very much reduced, since the ; 

largest portion of the flow can be determined from the simplified equations, and 
extremely fine resolution is only required for the breakdown region, where the full 
Navier-Stokes equations must be employed. 

The work carried out under this contract is concerned with the formulation and 
solution of the slender-vortex approximation and its numerical solution for 
incompressible and compressible flow; with the initiation of experimental work and 
the validation of the solution; and finally the prescription and determination of 
permissable inflow conditions for the slender-vortex approximation. The results of 
this analyses are presented in several publications (1] , (2) , (3] , 14] ; copies 
of the preprints and manuscripts are enclosed. : i 

Formulation of the Problem '> i 

i 

The slender-vortex approximation for compressible viscous flow was first formu- j 

lated in Refs. [1] and [2] . Therein it was shown that the flow is basicly < , 

governed by the local axial pressure gradient, which can be expressed through an >. • 

* 

integro-differential equation. According to this analysis the pressure gradient is j. 

mciniy determined by the radial distributions of the circumferential velocity j ! 

component, and the temperature. It could be shown that a hot core tends to • ' 

stimulate breakdown and that a cold core delays it. The direct influence of the ■ ’ 

; > 

compressibility is that the variable density tries to counteract the influence of . j 

positive axial pressure gradients, thereby delaying breakdown. . ; 



Permissnhle Inflow Conditions 


The slenderness condition imposed on the flow does not permit one to prescribe 
arbitrary inflow conditions, which are required for the integration of the simplified 
equations of motion. The inflow conditions are given by a set of radial profiles for 
the axial and circumferential velocity profiles, and the temperature. It w as shown 
in (11 and (2) that the inclination of the traces of the stream lines v/u in the 
axial plane satisfies a second-order nonlinear ordinary differential equation, to be 
integrated in the radial direction. The coefficients of that equation are solely 
determined by the inflow conditions. Since the inclination of the traces of the 
streamlines must obey the slenderness condition, the solution of the equation for 
v/u decides as to whether or not the prescribed inflow conditions for u, w, and T 
can be admitted. This is described in detail in Ref. [ 2 ) . 

Numerical Solution 


Finite-difference solutions of the slender vortex approximation were developed for 
incompressible flows, and compressible flows. The differential equations were 
casted into locally linearized implicit difference equations, which can be inverted 
by recursion. Details of the solution for incompressible flow are described in (4) ; 
the solution for compressible flow is described in [2] , and (3) . The most 
important results computed with the solution so far are also reported in the 
references just mentioned: In Ref. (4 ) , it is shown, for example, that the solution 
reacts sensibly to axial pressure gradients; comparison calculation were carried for 
the experimental data of Faler and Leibovich. The comparison, discussed in 
(4) shows that for vorticities with large swirl, the breakdown point cannot be 
predicted with the slender vortex approoximation, since it does not take into 
acocunt the upstream influence of the breakdown region. 

In Refs. 12) and 13) the results obtained so far for compressible flows are 
discussed. The influence of the freestream Mach number and of the temperature 
distribution prescribed for the inflow cross section on the flow development in the 
downstream direction is clearly demonstrated. 

Conclusions 


The slender-vortex approximation was analysed for incompressible and com- 




pressible flow. First the equations of motion were reduced in an order of magnitude 
analysis. Then compatiblity conditions were formulated for the inflow conditions. 
Thereafter finite-difference-solutions were constructed for incompressible and 
compressible flow. Finally it was shown that these solutions can be used to describe 
the flow in slender vortices. The analysis of the breakdown process must, however, 
be excluded, since its upstream influence cannot be predicted with the slender 
vortex approximation. The investigation of this problem is left for future work. 
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SCHLANKE WIRBEL IN KOMPRESSIBLER STRdMUNG 


E. Krause 

Aerodynamisches Institut der RWTH Aachen 
Wullnerstr. zw. 5 und 7, 5100 Aachen 


Der Einflu3 der Kompressibilitat auf schlanke Wirbel in stationnrer, reibungsbe- 
hafteter Stromung wird dargestellt. Eine entsprechende Ableittng fur reibungs- 
frcie Stromung ist in ( 1 ] und fur instationare Stromungen in l 2 J gegeben. 


Fur die Untersuchung wird angenommen, daG die Achse des Wirbels parallel zur 
ankommenden Stromung ist und daG die Stromung im Wirbel rotationssymme- 
trisch ist. Bezeichnet R den Kernradius in dem Querschnitt senkrecht zur 
Anstromrichtung, in dem der Wirbel erzeugt wird, und L die unbekannte Lange 
des Wirbels, gemessen in Stromungsrichtung, uber die der Wirbel nicht 
aufplatzt, so erfo rdert die Schiankheitsbedingung 

JL S 0( JjL)« i ‘ (1) 

Dabei ist v die rndiale und u die axiale Geschwindigkeitskomponente. Die mit 
Gleichung (1) vereinfachten Erhaltungsgleichungen fur Masse, Impuls und 
Energie entsprechen der Grenzschichtapproximation: 


Masse: -2— ( p u ) ♦ — — — ( p v r ) =0 

0 x r 0r 

(21 

Impuls, x-, r- und G -Richtung: 


0x 9r 0 x r 0r 0r 

(31 

pw^ _ 3p 
r " 0r 

(41 

»<■ i„if (4n 

9 x 9r r dr 0r r z 

(51 .. 

Energie 



Cplfu |I . »v |I I = „ |B. . v |H. . 1 JL ( r x|I| . , [[ r JL ,4 II 2 ] 

y 9x 9r 9x 3r r dr dr 9r 6r 


3r 


dx 


dr 


r dr 


( 6 ) 

Das System wird mit der Zustandsgleichung p = P/(RT) geschlossen, v/obei R 
die Gaskonstante ist. Zur Losung des Gleichungssystems (2)-(6) mussen Ein- 
strombedingungen 


= v„ . n<r : u< Mrl ws f»lr) . T = f-J r 1 


I7> 



I 


Symmetriebedingungen auf der Wirbelachse 

r = 0 ; x n < x : = v = w = = 0 , (8) 

u 3r 3r 

und Randbedingungen fur r 

r — .o , xq<x: usg^xl , w- 0 , T — g 2 (x) (9) 

vorgegeben werden. Die Funktionen f^(r) - f^(r) und g^(x) und g 2 (x) miissen als 
bekannt vorausgesetzt werden. 

Wegen der Schlankheitsbedingung (1) konnen f^(r) - f j(r) jedoch nicht willkurlich 
vorgegeben werden. Dies geht aus folgender Betrachtung hervnr; Aus der 
Kontinuitats-, Energie- und Zustandsgleichung laOt sich durch folgende 

Beziehung ausdrucken: 

SiL = . JL l£ . v H „ (Jidi w| v . _LJL| pvr )♦_!!_ (10} 

3x xp 3x T 3r r a 2 pr 3r pc p T 

Die Gro3e H steht zur Abkurzung fur die Dissipationsfunktion und den Warme- 
Ieitungsterm in der Energiegleichung. Nach Einsetzen von Gieichung (10) in die 
x-Impulsgleichung erhalt man nach Integration: 

• — = exp{-I)(/[(1 — ] exp ( 1} dr' 
u o o 2 pu 2 3x (11) 

-/I — (Di/pu - H/c D T)exp(I)dr*} 
o pu H 

Die GroGe 11^ steht zur Abkurzung der Schubspannung in der x-Impu.igleichung 
und I bedeutet 


i./iw-s (12) 

o a 2 r' 

Bei bekannten Einstrombedingungen f ^(r) - f ^(r) kann v/u nach Gieichung (10) 
bestimmt werden, wenn der ortliche Druckgradient 3P/3x unbekannt ist. Er 
laGt sich durch Differentation nach x und anschlieGender. Integration aus 
Gieichung (5) gewinnen: 


( x.r) = l x.r — oo ) 

3x ox 

♦/" VZr ) l A ( r'w ) - £ -Pr * 1 C — ) dr* 

c r 1 r‘ 3r T 3r* c p T r‘ u 

- 7 <3r*-/ => (2II- > — — H ) dr’ (13) 

r a 2 r' 3x r c p T r'u 


Die GroGe II 2 steht zur Abkurzung fur die Schubspannung in der 0-Impu!s- 
gieichung. 



Nacb Gleichung (13) ergibt sich stets ein nichtverschwindender Druckgradient 
— -•£ (x,r), da in reibungsbehafteter Stromung v/u stets von Null verschieden ist. 

O X 

Ein heiOer Kern kann dabei. den Druckgradienten vergroGern, wahrend die 
Umfangskomponente den Druckgradienten verkleinert. Dieser andert seinen 
EinfluG mit der axialen Machzahl. Das geht aus der x-lmpulsgleichung hervor, 
die sich in folgende Form bringen laGt: 


or u 



9p ,, w 2 v 1 1 9 . 3u , 

~ — -r) n — — - ) ♦ 

ox a* ur pu<- r or dr 


H 

pa*u 


(14) 


Gleichung (14) zeigt, daS bei Oberschalldurchstromung des Wirbels, d.h. u/a > 1 
der erste Term sein Vorzeichen andert. 


Sollen nur die Funktionen f^(r) - f j(r) auf ihre Kompatibiiitat mit der Schiank- 
heitsbedingung, Gleichung (1), im Einstromquerschnitt uberpruft werden, kann 
dies durch Elimination der Druckgradienten aus den Impulsgleichungen gesche- 
hen; Durch entsprechende Differentation nach x bzw. r und Subtraktion der 
resultierenden Ausdrucke voneinander, erhalt rnan nach Zusammenfassen 


( — J ♦ G« ( x ,r ) ( — 1 ♦ ( x.r ) f — )♦ G^(x,r ) = 0 

0r 2 u 3r u 1 u J 


(15) 


Die Funktionen G^(x,r) - Gj(x,r) lassen sich unmittelbar aus den Einstrombedin- 
gungen ermitteln. Integration von Gleichung (15) ergibt dann die radiale 
Verteilung v/u, so daG abgeschatzt werden kann, ob die Schlankheitsbedingung, 
Gleichung (1), erfullt ist. 
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NUMERICAL PREDICTION OF FLOW IN SLENDER VORTICES 


* • * 

Luis Reyna and Stefan Menne 

Aerodynamisches Institut, RWTH Aachen, West Germany 


Abstract 


We study the slender vortex approximation with attention put on high Reynolds 
number behaviour. It is shown that the breakdown of the approximation coincides 
with the criticality ' condition as introduced by Benjamin £ 12^ . We study free 
vortices with and without an adverse pressure gradient for viscous and inviscid 
flows. Finally we compare to experimental results from Faler and Leibovich . 


This research was conducted while the first co-author was at the Aero- 
\ 

dynamisches Institut as an Alexander von Humboldt-scholar. 




Correspondence and proofs for correction should be transmitted to Stefan 
Menne, Aerodynamisches Institut, Wijllnerstr. zw. Nr. 5/7, 5100 Aachen, West 
Germany. 
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NUMERICAL PREDICTION OF FLOW IN SLENDER VORTICES 



» ** 

Luis Reyna and Stefan Mcnne 

Aerodynamisches Institut, RWTH Aachen, West Germany 


1. Introduction 



Vortex breakdown was initially observed by Peckham and Atkinson ^ 1J for leading 
edge vortices formed on delta wings at large angle of attack and with large tip 
angles. The phenomenon has a drastic influence on the aerodynamical behaviour of 
the flow. In flows around wings its presence strongly decreases the lift £ and 
in combustion chambers it can be used to design flame holders £" 4 } . Despite of the 
large amount of research on this subject, the prob'em of vortex breakdown can not 
be yet considered as being fully understood. The present state of the art can be 
found in the review articles by Ludwieg [5] , Hall M and Leibovich £ 7^ . 

In experimental investigations Faler and Leibovich [” 8] found six different cases 
of breakdown for a swirl flow in a divergent tube. Two of the n are nearly 
axisymmetric and called bubble type, the remaining have either a spiral or heli- 
coida! shape. The spiral form is marked by a kink followed by a cork-screw shaped 
twisting of the vortex filament. In this case non-axisymmetric effects are im- 
portant and have to be included in the analysis of the flow £ . 



This research was conducted while the first co-author was at the Aero- 
dynamisches Institut as an Alexander von Humboldt-scholar. 

Correspondence and proofs for correction should be transmitted to Stefan 
Menne, Aerodynamisches Institut, Wijllnerstr. zw. Nr. 5/7, 5100 Aachen, West 
Germany. 



Vortex breakdown theories assume axial symmetry except Ludw.ey's £ lo] which 
regards it as a hydrodynamical instability of nearly axisymmetric flow. The 
important but unanswered question is the dependence of the breakdown type on the 
approaching flow. In order to explain the phenomenon Squire £ ll] , Benjamin 
[ 12,13 j and Bossel £ 14,15] introduced the concept of critical state. The flow is 
called supercritical when perturbations can only propagate downstream and subcritical 
when there is an upstream influence. Squire £ ll] suggested upstream perturba- 
tions would accumulate and produce breakdown at the critical state where the flow 
proceeds from super to subcritical and the phase velocity of the perturbations is 
zero. Benjamin £ 12] showed the existence of two equivalent solutions to a 
conjugate flow, one being supercritical and the other subcritical. He then regards 
breakdown as a sudden transition from super to subcritical states, similar to the 
hydraulic jump. It is not possible to predict vortex breakdown with this theory, but 
it allows a classification of flows. The significance of the theory can be seen in an 
experimental observation of Faler and Leibovich £ 8] : "All flows that exhibit 
vortex breakdown of the "axisymmetric" form (which we classify as types 0 and 1 
disturbance forms) or "spiral" form (our type 2) are supercritical upstream, in the 
sense of Benjamin." 

The details of breakdown can be studied by numerical integration of the full 
Navier-Stokes equations. Lavan, Nielsen and Fejer [ 16] , Kopecky and Torrance 
[17], Grabowski and Berger £ 18] computed axisymmetric, laminar, incom- 
pressible and stationary swirling flows. Grabowski and Berger £ 18] found a 
backward flow for subcritical initial profiles in contrast to the classification of 
Benjamin £ 12] . Since the double ring stucture was not present inside the bubble, 
Faler and Leibovich concluded due to their experiments £ 19] that the numerical 
solution should take into account time dependent periodic asymmmetric motions 
£19] . The time dependent calculation done by Shi £20] showed correctly this 
structure. The computed flow was stationary upstream of the breakdown point but 
iristationary and nearly time periodic downstream of it. This and similar 
computations are restricted to low Reynolds numbers, much lower than the ones 
present in technical applications. 

Flow pictures show nearly cylindrical stream surfaces upstream and this observa- 
tion can be used in order to derive an approximation of the Navier-Stokes equations 
usually called the slender vortex approximation. The assumption behind this 
approximation is analogous to the one from boundary layer theory. The slender 



vortex approximation has been used by Gartshore [21,22^ , Hali ^23-25j , Bosse! 
{^26^ , Mager [ 27j , Nakamura and Uchida [?&] and Shi £ 20^j . The 
approximation based on small gradients and small radial velocities fails at the 
breakdown point but is valid upstream of it and generally for stable vortices. 

The purpose of this paper is to study the high Reynolds number behaviour of the 
slender vortex approximation. In chapter 2 we prove that the breakdown of the 
approximation occurs at the critical state as remarked by Ludwieg 
chapter 3 we present numerical solutions for different Reynolds numbers. In the 
limit of no viscosity the external pressure gradient determines the breakdown. First 
preliminary results are shown in £ 29 J . Finally we compare to experimental data 
and discuss the advantages and limitations of the approximation. 


2. Slender vortex approximation 

Consider now the Navier-Stokes equations written in cylindrical coordinates 
(x, r, ©) with corresponding velocities (u,v,w) in a nondimensional form. The axial 
and radial velocity components are normalized with the axial velocity and the 
pressure with its value at the initial station for r-*co. The lengths are normalized 
with the vortex core radius which marks the region of viscous flow. The 
circumferential velocity is normalized with its value at the initial station at the 
edge of the vortex core. 

Including the slenderness condition 

v «“- » 

and assuming steady axial symmetric flow leads to a system of equations 


M.JX* + Uxcr -rp x f ( P'U-rV 

(Za.) 

T> — V / 1 

r r-~ 

(Zb') 

^ w* + tr £ (rw) r ** ~ r f ( r w) r ] r 

(2 c) 

\ 

+ (ru) r «■ O 

fcw) 


called the slender vortex approximation. Here the Reynolds number is based on the 
vortex core dimension and the undisturbed axial velocity and p is the pressure. 



I 


Notice that the type of this system is now/ parabolic for viscous flow and hyperbolic 
for inviscid flow coupled to two ordinary differential equations.’ The first such 
equation is the momentum equation in the radial direction and the second 

(nr) r f (r z u% ~ £-(?M r ) r ] (rxr) = 

, (3) 

can be easily derived from (2). 

Due to the slenderness condition the approximation based on small gradients and on 
a small radial velocity fails in the neighoourhood of the breakdown point. According 
to the theories of Squire Cn] and Benjamin the flow is supercritical 

upstream of this position and perturbations can only propagate downstream. As a 
consequence the influence of the breakdown bubble is not contained in the slender 
vortex approximation. 


Symmetry conditons are imposed along the r = 0 axis: 

AjLf- o , \r~o , w ~ O faa) 

At the outer radius the type of the system allows three boundary conditions for 
viscous flow and one for the pressure for inviscid flow; the physical boundary 
conditions for free vortices are 

(?+i*M z )(rx)-^ c t VJ—+0 ? — (tb) 

V-> CD r - > “ *—+CD 

where C is a constant taken from the initial conditions and p. (x) is the external 
pressure. These boundary conditions are also valid for the inviscid system when the 
radial velocity has a negative sign at infinity, otherwise only the pressure 
can be given 


f—y co 


Otc) 


The slender vortex approximation in vorticity, circulation, stream-function 
formulation becomes 


a. +■ r ^( r ) r “ r 2, T x ^ ( r 

mV* + vT r (£r r ) r 
-a=(fTr) r 


(So.) 


(Sb) 

(Sc) 



with the vorticity /2 = -u r , local circulation T 7 = rw, stream-function ^ nr.-! 
velocity compor „-nt.s 

IM 

(it) 

In the vorticity, circulation, stream-function formulation the swirl influence is 

i "P 2 

given explicitly in the vorticity transport eq''"*ion by the term / 3 I x in contrast 
to the formulation in primitive variables where this coupling follows implicitly ov.r 
the pressure. For an isolated vortex it is advantageous to consider the perturbation 
of the stream-function from parallel flow. In this case the parameter £ in eq.(ba) 
is equal to one; else 5 ir equal to zero. 

The boundary conditions (4) corresponding to the vorticity-stream-function for- 
mulation are at the axis of symmetry 

12=0 , V = 0 , O (7a) 

At the outer radius t'..e boundary conditions (4b) translate into 

.a — , V — >(!>>% — > r(4-h(c-? co (x'})') W 

r— >oo r ■ — > do r — v co 

1 2 

where C =.(P. + ^ u. ) is a constant taken from the initial conditions and (x) is 
the external pressure. The parameter specifies the rate of the circumferential 
velocity at the edge of the vortex core to the axial velocity at infinity in the initial 
station. T hese boundary conditions are also valid for the inviscid system when the 
radial velocity has a negative sign at infinity, otherwise only the pressure 
can be given 

y — (?c) 

r-*co 

Finally for viscous flows in a pipe the no-slip condition 

, TT = O , W - O 

holds at the outer radius. 

For free vortices the corresponding boundary conditions for (3) arc 

V=0 ext r=0 CLKcL V r t = 0 ^ r — >CD (?0 

where u^Cx) = ZfC-P^, (x)) is a given function. 





Initial conditions have to be provided for the axial and circumferential velocity 

'“‘(r, x a ) - <r> 

ftf) 

= Wo cn 

The pressure, the radial velocity, the vorticity, the circulation and the stream- 
function can be obtained from them. 

Solutions of the slender vortex approximation con cease to exist due to two 
different mechanisms: first, when the axial velocity vanishes somewhere in the 
flow and second, when the radial velocity becomes unbounded. The first case is 
similar to the boundary layer separation while the second is related to the 
criticality condition introduced by Benjamin [ 12 J . The second situation arises 
since there is no viscosity present along the axial direction and therefore 
unbounded gradients along it can appear. Since the viscosity present in (2) only 
controls the radial gradients of the axial and circumferential velocity components 
axial gradients and also the radial velocity can grow infinitely. Moreover, in this 
case the solution tends to that of the homogeneous problem (3). 


(r V> rr -f(rv) r r f-^rr (r V) r -£ (f «r) r ] <r v) = O 



Indeed numerical solutions of the slender vortex approximation shnw that the 
viscous forces are really small compared to pressure and inertia forces near 
breakdown of the slender vortex equations. 

As in boundary layer theory the behaviour of the approximated and the full Navier- 
Stokes system are identical shortly before breakdown when the proper boundary 
conditions are imposed. We return to the point v/hen comparing predictions of the 
slender vortex approximation and experimental results. 

We now compare this breakdown - condition to the criticality condition as introduced 
by Benjamin £l2^ (see also Hall [6^] )• 


Benjamin considers steady, inviscid and axisymmetric flow. Introducing the 

' 1 1 12 2 2 
streamfunction if/ with u = - and v = - -yv, the total head h = p + ^{u +v +w ) 

and the circulation "P = rw the equation of motion for these flows can be 

summarized in a single equation for the stream-function. 


dx r 


ipi = . r ir. tr ^ 

r sr 1 dy dy 


w 


do) 


We now consider a perturbation of a quasicyiindrical flow 

Y(x, i') = Y(x,r) + £ ^(x,r) e rx 

where £« 1 and “ 9 “^ ^ "g~- Th 's perturbation introduced in ( 9 ) and neglecting 
higher order terms in £ yields an equation for F. The critical state F corresponds 
to y = 0 and in this case 

Tc Fr -fH r +r-£(?A r ) r + 7^T(rV) r ]F c =0 (11) 

If the solution with initial data = 0 and F cr = 1 at r=0 vanishes away from the 
axis the flow is subcritical; when this solution vanishes at both the axis and outer 
radius R it is said to be critical and otherwise supercritical. 

In comparison to eq.(ll) the equation for the radial velocity ( 8 ) 

(rv) r + r-£(f« r ) r t T ^ r -(r*w») r ](rv) = 0 (- 11 ) 

shows an identical form with the boundary conditons 
V*0 at r«=0 and tT r + ~ + - c ^ op - = O at fa) 

where the outer radius R has a finite value with R»l. 

The general solution of eq.(12) reads 

(rv) = C H(r) (14) 

where the boundary condition ( rv )( r _gj = 0 is already used. 

With the constant C the solution function can be adapted to the boundary 
condition at r=R. Excluding the trivial solution vsO, i.e. v^_ R) £0, the demand 

r -* it •• h — » 0 

leads to a critical case, since necessarily the constant C has to go to infinity and 
consequently 

o< r< n •• (rv) — > 00 

This characteristical behaviour of the radial velocity is present in the slender 
\ 

vortex approximation in the vicinity of the breakdov.m region. A comparison of this 
critical case with the previously described critical flow state according to the 
theory of Benjamin [ 12*] leads to a complete agreement: eq.(ll) and (12) are 
equivalent to each other just as the boundary and the criticality conditions since 
the condition F cr = 1 excludes only the trivial solution. 




Now we con see that the breakdown condition for the slender vortex approximation 
corresponds to the condition from Benjamin [ 12^ and breakdown can be seen as a 
transition from supercritical to suberitical flows as remarked by Ludwieg [ 5^ . 
Hall £ 6]] reached also the same conclusion. He explained the equivalence using 
the existence of two solutions for the slender vortex a small distance apart from 
each other which do not coincide os the distance tends to zero. He then shows that 
"the condition for the appearcnce of arbitrarily large axial gradients turn out to be 
identical to the condition for the critical" [6]. 


3. Solution and results 

We used both the primitive variables and the vorticity-stream function formulation 
for the numerical solution of the slender vortex approximation. In a first approach 
we used primitive variables with centered discretization. The radial velocity is 
evaluated only at the midpoints of axial intervals (see fig. 1). The discretization 
only use two "time" levels making it convenient for variable axial spacing. The 
function value distribution is shown in fig. 1. 


Here fj. denotes the numerical solution at the (i,j) node corresponding to 
(x,r) = (iAx,jAr), where Ax and A r are the axial and radial spacing. The equations 
are solved marching along the axial direction as the type of the system indicates it. 
At each new axial station a system of nonlinear equations must be solved. This is 
done using Newton's approach and a linear band solver for the Jacobian inversion. 


In order to decrease the band width, the symmetry condition for the axial velocity 
is applied as 


/-Ly ~ 


AT 


AT 


Ke[ rJTv? .1 -v q (ar 1 ) (IS-) 

1 ZAX &X J ' 


l AX 




instead of three point extrapolation formula. (This boundary condition also holds in 
the inviscid cose). The discretization is second order accurate in both radial and 
axial direction and unconditionally stable with respect to the size of Ax. 

Nevertheless during the calculation the size of /^x is decreased when the Courant- 

i A y 

Friedrichs-Lewy number CFL = max. — exceeds a predetermined constant 

CFL . A x is also decreased when tfie Newton's procedure fails to reduce the 
max 

residual below the tolerance limit TOL after NEW iterations. 

max 

-A 

Normally we use CFL =2, NEW =3 and TOL=10 , being the solution not 
max ’ max 5 

sensitive to these values. 


The second numerical approach was applied to the vorticity stream-function 
formulation (6) . We use again centered discretization with function values 
determined at node points except the radial velocity which is evaluated at inter- 
mediate axial points (i+ ^,j) (Crank-Nicolson formulation). The resulting discreti- 
zation is unconditionally stable and second order accurate. At each new axial 
station a tridiagonal system must be solved for the stream-function, vorticity and 
circulation. For a free vortex the equations are not coupled through the boundary 
conditions and therefore can be solved iteratively as follows. We initialize 

(o') 

<<» J m 

and using these values compute T ^ and theni^f^* 

Next the stream-function^^ is computed and a new iteration can be started. The 
procedure is repeated until 
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where £ is usually set to l(f ^ 10 . The iterative procedure has a low storage 
requirement, only two axial stations has to be kept. It has the drawback to be only 
convergent for, supercritical profiles. This is not the case for the Newton’s-iteration 
approach that allows subcritical initial profiles. On the other hand initially 
subcritical profiles are of little physical interest since these flows con be regarded 
as already broken down according to the theories of Squire [ llj and Benjamin 

M: 








I 


The convergence problem translates into increasing number of iterations as the 
flow reaches criticality conditions. In these situations the axial step is decreased in 
order to keep the amount of work low; this also gives a better axial location for the 
breakdown point. Each line the axial step is decreased, the iteration is restarted 
with the last converged values. 


Since few grid points are required outside the vortex core we use a transformation 
in the radial direction 

r : — * or 

_r_ = <an (J g? 

The uniform spacing in the £ -direction gives almost a uniform distribution in the 
radial direction when G’ ma „«l and an accumulation of grid points near the axis 
when GT max ~ 1* The pressure is obtained integrating the radial momentum 
equation 

?i-j - T+.j« 2i-< + + if <9 o + )] ^ 


where 
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The initial value is the pressure at O' . g is an even function of 6 . In order to 

max 3 

decrease the truncation error the pressure is integrated analytically outside the 
vorticity core where the circulation is constant. 


The initial values are taken as 
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These polynomial distributions are already used by Mager £.27} , Grabowski and 
Berger £l8|] 'and Shi £.20^ snd "...were chosen to approximate the experimen- 
tally-measured velocity in vortex cores such as those of trailing vortices... " [18] • 

The parameter c<, controls the shape of the axial profile, uniform flow for 0C = 0, 
jctlikc foroC>0 and wakelike for 0. The circulation corresponds to solid body 
rotation for r« 1 and to potential flow outside the core of vorticity. The 
polynomials fulfill symmetry conditions and provide a smooth transition form the 



solidly rotating core into the outer flowfield. The parameter (3> specifies the rate 
of the circumferential velocity at the-edge of the vortex core to the axial velocity 
at infinty in the initial station. 

We now return to the boundary conditions for free vortices. We assume constant 
pressure gradient 

£Es£ s 7 -* (Zi) 

dx 0 

Then the slender vortex equations can be scaled using 

X-X't , - ,• 

and therefore the pressure gradient for the transformed system is equal to one. The 

1 

flow is relevant for values x < ■= since the outer continuitly desacceleration of the 

^ ~ 1 
outer flow makes the outer axial velocity vanish at x = ^. 


X = 


(23) 


When no pressure gradient is present in the outer flow the scaling 
X V 

makes the equations and the initial conditions Reynolds number independent. 
Notice that the axial coordinate is stretched but the radial coordinate is left 
unchanged in contrast to boundary layer type of scaling (e.g. Hall £,253 )• The 
result is an increasing breakdown length for increasing Reynolds number in contrast 
to experimental observations £303 for which the location reaches a limit. Due to 
this behaviour at high Reynolds number we do not expect this to be the physical 
mechanism behind breakdown which we believe to be pressure induced. 


For tube flows the numerical procedure for the stream-function formulation was 
slightly changed due to coupling of the vorticity and stream-function through 
boundary conditions. The stream-function was computed from the fourth order 
differential equation obtained when (6c) is introduced in (6a) and using a five point 
centered discretization. For the primitive variable formulation these new boundary 
conditions introduced no new complications. 

In order to resolve the boundary layer on the tube wall, the following stretching 
function was used 

r = 'R'(x) Oyrcia.y\. (<3 (an. & 4) 

with being the tube radius and = 0.0. 



Paler and Lcibovich [*] measured several velocity distributions which were used 
as initial condition for the flow calculations. 


We study our breakdown criterion described in chapter 2 by applying it to free 
flows driven by a pressure gradient which reach critical state already in the initial 
section. We used simplified profiles that allow us to do analytical work: 
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These profiles ore an approximation to the profiles used in the numerical examples. 
The axial velocity distribution is discontinous and the circumferential velocity can 
be described by solid body rotation inside the core and potential flow outside of it. 


Introducing these distributions in (12) we obtain 
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Since the axial velocity is discontinous jump conditions have to be prescribed 
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where CO denotes the jump in the values of f across r. 
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The solution of eq.(16) reads 

v = c A y* (z) 

or = C z r + yr 

where is the Bessel function of the first kind of first order. The corstant c^ is 
determined by the pressure gradient at infinity. Since the condition v = 0 at r = 0 is 
already used, the constants c^ and c^ are determined by the jump conditions (27). 

It follows 
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Criticality is already reached when 

?o |V?) = o (30) 

that is 

<x & = -T^—P>-'i witk j.0 K i.wr <»3 

^ 3°'° ® 

(j is the first non trivia! zero of the zeroth Bessel function 3 ). This critical 
curve is drawn in figure 2 as a dotted line. The curve marked by triangles shows the 
numerical profiles which are initially critical. Due to the undershoot in the axial 
velocity for c£< 0 of the discontinuous analytical distribution when compared to the 
numerically used profiles, the limiting curve lies too high producing a destabilizing 
effect (tf>0 produces the inverse effect). The overshoot of the analytically used 
profiles for the circulation shifts the curve &(-.((*>) upwards. The numerical results 
can be summarized as follows: Or the right side of the limiting curve in fig. 2 no 
solution of the slender vortex approximation can be found (subcritical region). At 
the left there exists a region where vortices have a limited breakdown length. In 
the following region at the left no vortex breakdown appears at all. The vortices 
are stable and dissipate themselves away. 


3.1. Isolated vortex for inviscid and viscous flow 

Figs. 3-6 shovy results for the free vortex with CC = 0, (i= .8 and ^ _ 0 . This 

' x r=oo 

vortex breaks down at ^ = 0.013. In fig. 3 we show the relative strength of the 

forces present in the axial momentum balance. Here Fp is the pressure force, Fj 

the inertia and Fy the viscous force. Initially pressure and inertia forces are almost 

equal while viscous forces are much smaller. Near breakdown both pressure and 

inertia forces show a dramatic inc'easc in their strength. This increase is not 




present in the viscous forces. Near breakdown the flow is essentially inviscid, 
viscosity plays a secondary role. This is in agreement with the corresponding 
theories of Squire [ll"J and Benjamin £l2^ and the critical condition which calls 
for a non trivial solution to the homogeneous version of equation (3). 

Fig. h shows the radial dependence of the three components of the velocity vector 
at different axial stations. The circumferential velocity decreases in magnitude and 
• the axial velocity develops a wakelike profile as the flow reaches breakdown 
conditions. The radial velocity shows a dramatic increase just before breakdown, 
similar to that of the pressure and inertia forces. Fig. 5 presents the radial distribu- 
tion of the pressure. The behaviour of the pressure is quite similar to that of the 
axial and circumferential velocity component. Fig. 6 shows the axial distribution of 
the pressure and of the radial velocity at different . adial distances from the sym- 
metry axis. Notice an almost linear change in the pressure until shortly before 
breakdown. In fig. 7 and 8 the velocity components and the pressure are shown at 
different axial stations for inviscid flow. Fig. 9 demonstrates the axial dependence 
of the pressure for inviscid flow. The good agreement of the breakdown process 
compared to the viscous case indicates that breakdown is an inviscid phenomenon. 

Fig. 10 shows the dependence of the breakdown length on the initial profiles for 
free vortices without adverse pressure gradient. Here only supercritical profiles are 
considered. We also show the corresponding results obtained by Shi [ 20^ by 
solving the time dependent Navier-Stokes equations. The slender vortex approxima- 
tion is in good agreement with his numerical results. 

For a given [b, there exists a lower bound for &■ leading to initially supercritical 
profiles. It is possible nevertheless to numerically solve system (2) starting with 
subcritical profiles when using primitive variables in conjunction with a direct 
solver. V/e found that all such flows exhibit breakdown at some axial distance from 
the initial profile. 

The increase of the jet type profile provides a stabilizing effect on the flow. For & 
large enough and passed a critical value the vortex does not break down at all and 
viscosity only flattens the profiles (this is usually referred as aging of the vortex) 
(see fig. 11). The effects of 0£ on the flow are still important near its critical 
value. For <X passed this value the axial velocity at the axis decreases initially but 
after reaching a minimum value increases again with no breakdown of the flow. 

There is no unbounded growth in cither the pressure gradient or the radial velocity 
any longer. 





















Fig. 12 shows results for vortices with constnnt external pressure gradient 

^ ^ and fig. 12a an increase in the breakdown length for increasing 

r= oo 

Reynolds number with its value attaining a limit for Re-»co. This is the expected 
behaviour for vortices that exhibit breakdown without external pressure gradient. 
Vortices that are stable in absence of pressure gradient (here oC = 1 . 1 ) can break 
down for Re 'X large enough (in this case Re*T^ 0.1) but are stable otherwise, as 
once more expected. 

Fig. 12b shows the dependence of the breakdown length for small pressure gradient. 
For T small we can see x^^/Re reading the expected values obtained from the 
vortex without pressure gradient shown in fig. 10. For increasing Reynolds numbers 
the ratio x^p/Re tends to zero since x^p tends to the limit shown in the left 
picture. Fig. 13a and 13b show similar results obtained from different swirls. We 
can see again the disappearence of breakdown for stable vortices for small enough 
pressure gradients. 

In order to stress the influence of the pressure gradient on the flow we sketch in 
fig. 14 the results from figs. 10, 12 and 13 combined for Re *T = 10. The flow is 
shown to be extremely sensitive to small pressure fluctuations at large Reynolds 
numbers. A change in 1% in the pressure, normalized with the dynamic head, over 
one vortex core length translates into a 100% change in the breakdown length for a 
Reynolds number of 1000. 

The behaviour of the flow shown in figs. 12 and 13 can also be found in 
experimental investigations. Werl6 [[ found that the breakdown length of a 
leading edge vortex formed on a delta wing does not change for high Reynolds 
numbers (Re^> 10**) for an angle of attack of 20 degrees. If we assume that the 
adverse pressure gradient of the vortex lies approximately in the range of 0.5 to 2% 
of the dynamic head over the vortex core radius, then the value of RfT lies 
between 50 and 200. 

For these values the breakdown length has approximately reached its ultimate 
value accepted for Re-T— » co (Fig. 12a and 12b). The extreme sensitivity against 
adverse pressure gradients shown in fig. 14 can be observed also in fig. 15 (from 
[30] ). Although the obstacle is positioned f3r downstream of the breakdown 
point, there is a significant change in the breakdown length. 
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3.2. Comparison to experimental results 


Faler and Leibovich O] performed several experiments on vortices in divergent 
tubes. The swirl component was generated by swirl vines and afterwards the flow 
was led through the tube. Profiles for the axial and circumferential velocities were 
measured diameter upstream of the divergent part of the tube. The opening 
angle for the tube was 1.43° and various mass fluxes and swirl vane angles of 
attack investigated. They found six different breakdown types for the flow, the 
corresponding breakdown lengths and their dependence on the Reynolds number and 
swirl parameter are shown in fig. 16. The Reynolds number is based here on the 
averaged axial velocity and the diameter of the tube. 

Since the slender vortex approximation only deals with axisymmetric flows, only 
bubble-type breakdown is considered in this comparison (types "0" and "1" according 
to tne terminology used in £ b"] ). Stable bubble structures are only present for 
certain ranges of circulation 61 and Reynolds number Re. 

One more inconvenience causes from the lack of detailed information on the 
boundary layer structure for high Reynolds numbers. Fig. 16 shows the numerical 
and computational results combined. The initial distributions are constant along 
each solution curve and coincide with' the experimental values only at nodes marked 
(*). For (O = 1.07 we can see a good agreement, the trend of decreasing breakdown 
length for increasing Reynolds numbers is present and the numerical values are 
similar. For Re around 2300 and below no breakdown occur in the numerical 
experiments. 

For 60 = 1*56 the . .merical prediction of the breadown length decreases for 
increasing Reynolds number larger than 3000. This behaviour reverses itself for 
1000 5 Re £ 3000. For Re ^ 1000 there is no breakdown present in the numerical 
experiments. The prediction of the breakdown length is now far off from its 
experimental counterpart. The physical reason behind this phenomenon is the 
sensitive behaviour of the flow on pressure variations. The pressure prediction 
obtained from the slender vortex approximation does not take into account the 
influence of the bubble presence on the pressure field acting on the vortex core. We 
propose that a good prediction of the breakdown location can be obtained from an 
inviscid calculation for the vortex core with an externa! pressure gradient obtained 
from the numerical solution of the full Navier-Stokcs equations without the 
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slenderness approximation. This is an explanation for the fact that the numerical 
prediction strongly deviates from the experimental results. Moreover, the assump- 
tion of selfsimilar inflow profiles is not fulfilled in this case (Re = 6000). To 
compare with experimental data we had to assume selfsimilarity to transfer to 
other Reynolds numbers because only a few inflow profiles are given in [bJ . 

t 

The lack of upsteam influence in the approximation leads to predictions for the 
breakdown length which are longer than the observed values. Nevertheless, for all 
calculations of (0 and Re the inflow profiles are supercritical and approach critical 
state near breakdown. 

Faler and Leibovich make the following remark: 

"All flows that exhibit vortex breakdown of the "axisymmetric" form (which we 
classify os types 0 and 1 disturbance forms) or "spiral" form (our type 2) are 
supercritical upstream, in the sense of Benjamin" £8^ . 


A. Conclusions 

The slender vortex approximation was studied in particular for high Reynolds 
numbers. For free vortices without external pressure influences the breakdown 
length is proportional to the Reynolds number. For free vortices with adverse 
pressure gradients, the breakdown length is inversely proportional to the value of 
its gradient. Flows with small pressure gradients take a long distance t" breakdown. 

For low Reynolds number the prediction of the simplified system agree quite well 
with the. ones obtained from solutions of the full Navier-Stokes equations. It was 
found that the flow becomes quite sensitive on pressure fluctuations for high 
Reynolds numbers and that the failing of the slender vortex equations corresponds 
to the critical condition from Benjamin [] 1 2^j for inviscid flows. The last comment 
holds since viscous forces are negligible near breakdown compared to inertia and 
pressure forces. The viscous forces do play a role for low pressure gradients by 
controlling the aging process. The predictions from the approximating system were 
compared to experimental results and for low swirl a good agreement was obtained. 
For higher swirl upstream effects on the pressure produced by the breakdown 
bubble deteriorates their agreement. This can not be incorporated into the slender 
vortex approximation. 
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Tig. 1: 


Coordinate system and grid definition 


Fig. 2: 


Stability diagram for vortex breakdown 

Transition from super- to subcritical flow according to Benjamin 
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Breakdown at inflow section for slender vortex approximation, initial 
condition (20) 

Region of temporary stable vortices 

Breakdown at inflow section for slender vortex approximation, initial 
condition (25) 


Fig. 3: Axial variation of the inertia, pressure, and viscous forces per unit 

volume as computed with the viscous slender-vortex approximation 


Fig. h: Radial profiles of the axial, radial, and circumferential velocity 

components computed with viscous slender-vortex approximation 


Fig. 5; Radial profiles of the static pressure as computed with viscous slender- 
vortex approximation 

Fig. 6: Axial variation of the axial velocity component and of the static pressure 

as computed with the viscous slender-vortex approximation 


Fig. 7a: Axial velocity profiles computed with slender vortex approximation for 

inviscid flow 


Fig. 7b: Radial velocity profiles computed with slender vortex approximation for 

inviscid flow 


Fig. 7c: Circumferential velocity profiles computed with slender vortex approxi- 



Fig. 0: 


Fig. 9: 


Fig. 10: 


Fig. 11: 


Fig. 12a: 


Fig. 12b: 


Fig. 13a: 


Fig. 13b: 


Fig, 14: 


Pressure profiles computed with slender vortex approximation for 
inviscid flow 


Axial pressure variation computed with slender vortex approximation for 
inviscid f'nw 


Shape parameter Ci. as a function of the breakdown length. Computed for 
several swirl rates 

Axial variation of the axial velocity component computed with viscous 
slender vortex approximation 


Breakdown length as a function of the Reynolds number for constant 
externally imposed pressure gradient T. Computed with viscous slender- 
vortex approxiamtion 

Breakdown length as a function of. imposed external pressure gradient X 
for constant Reynolds number. Computed with viscous slender vortex 
approximation 


Same as Fig. 12a. The swirl parameter (b - 0.8 


Same as Fig. 12a. The swirl parameter (b = 0.8944 
\ 

Influence of an externally imposed axial pressure gradient X on the 
breakdown length os computed with the viscous slender-vortex approxi- 
mation 
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Fig. 15: Influence of an externally imposed axial pressure gradient on the 

breakdown length. (From £ m3 ) 


Fig. 16: Comparison of experimental data (from £ 8^ ) and numerical results with 

slender vortex approximation. The type of disturbance (0-6) and its mean 
axial location vs. Reynolds number. 






